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We apply a complex chiral random matrix model as an effective model to QCD 
with a small chemical potential at zero temperature. In our model the correlation 
functions of complex eigenvalues can be determined analytically in two different 
limits, at weak and strong non-Hermiticity. We compare them to the distribution 
of the smallest Dirac operator eigenvalues from quenched QCD lattice data for 
small values of the chemical potential, appropriately rescaled with the volume. 
This confirms the existence of two different scaling regimes from lattice data. 



1. Introduction 



The presence of a chemical potential /i in QCD remains a challenging prob- 
lem, due to the fact that the Dirac operator eigenvalues become complex. 
One possibility to study this situation analytically is given by random ma- 
trix models. The origin of such an effective description for real Dirac eigen- 
values at fi = 0, as first suggested in [1], is by now very well established, 
and we refer to [2] for a review. A chiral random matrix model including 
a chemical potential was proposed in [3]. While its partition function and 
susceptibilities were determined analytically in [4], the microscopic corre- 
lation functions are not known to date. Very recently the model has been 
simulated on the lattice [5], confirming the results of [4]. In [6] complex 
Dirac spectra have been analyzed using quenched lattice data in the bulk 
of the spectrum. Here chiral symmetry is not important, and the relevant 
matrix model correlations are given by the known results for the Ginibre 
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ensemble [7]. A transition starting from the Gaussian unitary ensemble 
(GUE) at fi = via the Ginibre to the Poisson ensemble was observed [6]. 

Very recently an alternative chiral matrix model with complex eigen- 
values was introduced and solved in [8]. Here the microscopic, complex 
correlation functions are calculated in two different limits. The first limit is 
called weakly non-Hermitian and was discovered in [9]. It interpolates be- 
tween the chiral GUE [1] and the correlations in the second limit at strong 
non-Hcrmiticity. The results of Ginibre [7] are at strong non-Hermiticity, 
too, but differ from [8] due to chiral symmetry. We note that the model [8] 
is always in the phase with broken chiral symmetry, in contrast to [3] . 

2. Chiral random matrix theory with complex eigenvalues 

In this section we summarize the results of the complex extension [8] of the 
chiral GUE. The matrix model partition function is defined as 

N 

Z^\t) = / \{dz 2 |z,| 2a+1 exp 

*ni*2-*?r- w 

k>l 

In analogy to the chiral GUE [1], the power a — Nf + v combines the 
number of massless quark flavors Nf and the topological charge v. The 
non-Hermiticity parameter r G [0, 1] relates to the chemical potential \x via 

V = 1-t 2 (2) 

from comparing [3] and [8] at small fi. For r — > 1, the complex eigenvalues 
Zj become real, and we recover the partition function of the chiral GUE [1]. 
In the limit r — > 0, the non-Hermiticity is maximal, and eq. (1) becomes 
a chiral extension of the Ginibre ensemble [7]. The relation between the 
matrix model [3] and eq. (1) consists in a nontrivial mapping of the complex 
phase of the Dirac determinant in the former to the exponent of the latter, 
and we conjecture that the two models are in the same universality class. 

The model of eq. (1) is solved at finite and large N using the technique 
of orthogonal polynomials in the complex plane [8]. We distinguish two 
different large- TV limits: in the weakly non-Hermitian limit [9] we keep 

lim lim N(l - t 2 ) = a 2 (3) 

jV— >00 T— >1 

fixed. We thus send fj, — > such that the volume V times [i 2 is fixed. 
In this limit the eigenvalues macroscopically collapse to the real line while 
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microscopically they still extend into the complex plane. At strong non- 
Hermiticity we take N — > oo at fixed r € [0, 1). In this limit the macroscopic 
eigenvalue density becomes constant on an ellipse in the complex plane. 
The result for the microscopic density at weak non-Hermiticity reads a 

/ dt e- a2t J a (VtOJ a (VtC) , (4) 
Jo 

where we have rescaled the complex eigenvalues £ = Nzy2 with the same 
power in N as for real eigenvalues [1]. In contrast, we obtain for the micro- 
scopic density at strong non-Hermiticity a 

Pstrong(0 = jj |f| exp 

rescaling £ = y/2N z here. For the higher order fc-point eigenvalue correla- 
tion functions in both limits we refer to [8] . 



erf(a) 



Ifl exp 



-1 
2(1 -r2) 



la 



2(1 -r 2 ) 



(5) 



3. Comparison with lattice data 

In this section we compare quenched QCD lattice data with the predictions 
of eqs. (4) and (5). The data were generated on a 6 4 lattice at gauge 
coupling (3 = 5.0 using staggered fermions. We have chosen two different 
values for the chemical potential, fi — 0.006 and /i = 0.2, with 3500 and 
4500 configurations, respectively, corresponding to the two different regimes 
of weak and strong non-Hermiticity. In both cases, we are concentrating on 
the microscopic scale (i.e. a few eigenvalues) so that instead of unfolding 
the data we only need to rescale them by the mean level spacing. 



1.4 
1.2 : 
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Figure 1. Cuts of the Dirac eigenvalue density for fixed real part of the eigenvalues, plot- 
ted vs the imaginary part, for /i = 0.006, along with eq. (4). Left: real part corresponds 
to first maximum in Fig. 2, right: real part corresponds to first minimum in Fig. 2. 



a Note the different normalization compared to [8]. 
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Figure 2. Left: cut of the Dirac eigenvalue density along the real axis for fi = 0.006, 
along with eq. (4). Right: the density in the complex plane, along with eq. (4) (which is 
slightly displaced for comparison). 

For fi = 0.006 we determine the parameter a from the weak limit by a 
fit to the decay into the complex plane for a fixed real value, as shown in 
Fig. 1. We find on average a = 0.19, varying about 4% when determined 
along the first maximum or minimum on the real axis. This agrees roughly 
with a ~ 0.27 obtained from eqs. (2, 3) and the rescaled level spacing. 
The discrepancy may be due to the uncertainty in the latter and will be 
further investigated. In Fig. 2 (left) we show a section of the data along 
the real axis. Since a is small, we could equally well plot the real density 
Mi = l^iMO 2 + MO 2 }, differing by < 0.5%. The histogram in 

Fig. 2 (right) shows the quantitative agreement in the complex plane. 

At strong non-Hermiticity we rescale the data for (i = 0.2 independently 
in the real and imaginary directions with the square root of the respective 
level spacing. Cuts along the axes are shown in Fig. 3, and the complete 
plot is in Fig 4. Here we have used the relation (2) instead of fitting, leading 
to t = 0.9165. The repulsion of the eigenvalues from the origin, which is 




-1 -0.75 -0.5 - 



Figure 3. Cuts of the Dirac eigenvalue density for /i = 0.2, along with eq. (5). Left: 
cut along the imaginary axis, right: cut along the real axis. 
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Figure 4. Dirac eigenvalue density and eq. (5) in the complex plane for fi = 0.2. 

very different from that in the weak limit, is clearly seen in the data. 
4. Conclusions 

In this first exploratory study, we found good agreement between the pre- 
dictions of the model (1) and data from lattice QCD simulations. In partic- 
ular, the two predicted scaling regimes corresponding to weak and strong 
non-Hermiticity are clearly visible in the data. Obviously, more work is 
needed to make these conclusions more quantitative (e.g. larger volume, 
higher statistics, unfolding details). 
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